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ABSTRACT 

We analyse the far infrared dust emission from the Galactic Centre region, including the 
Circumnuclear Disk (CND) and other structures, using Herschel PACS and SPIRE photometric 
observations. These Herschel data are complemented by unpublished observations by the In- 
frared Space Observatory Long Wavelength Spectrometer (ISO LWS), which used parallel mode 
scans to obtain photometric images of the region with a larger beam than Herschel but with a 
complementary wavelength coverage and more frequent sampling with ten detectors observing 
at ten different wavelengths in the range from 46 to 180 /im, where the emission peaks. We 
also include data from the Midcourse Space Experiment (MSX) at 21.3 /mi for completeness. 
We model the combined ISO LWS continuum plus Herschel PACS and SPIRE photometric data 
toward the central 2 pc in Sagittarius A* (SgrA*), a region that includes the CND. We find that 
the FIR spectral energy distribution is best represented by a continuum that is the sum of three 
grey-body curves from dust at temperatures of 90, 44.5, and 23 K. We obtain temperature and 
molecular hydrogen column density maps of the region. We estimate the mass of the inner part of 
the CND to be ~ 5.0 x 10 4 M Q , with luminosities: L cavity ~ 2.2 x 10 6 L© and L C to ~ 1-5 x 10 6 
Lq in the central 2 pc radius around SgrA*. We find from the Herschel and ISO data that the 
cold component of the dust dominates the total dust mass, with a contribution of ~ 3.2 x 10 4 
M Q ; this important cold material had escaped the notice of earlier studies that relied on shorter 
wavelength observations. The hotter component disagrees with some earlier estimates, but is 
consistent with measured gas temperatures and with models that imply shock heating or turbu- 
lent effects are at work. We find that the dust grain sizes apparently change widely across the 
region, perhaps in response to the temperature variations, and we map that distribution. 

Subject headings: -dust, extinction -Galaxy: center -infrared: ISM -ISM: individual objects (SgrA*) 
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1. Introduction 

The central — 1 parsec of the galaxy is a low 
density central cavity that surrounds SgrA*, the 
compact variable radio source located close to the 
dynamical centre of the Galaxy - a black hole 
of ~4.5xl0 6 M^ at a distance of 8.3±0.4 Kpc 
dKerr et alJll986t iGhez et alJl2008t iGillesen et al ' 
20091 ) . The central cavity is characterized by emis- 



sion of low ionisation atomic species. The cavity 
contains a cluster of Wolf-Rayet (WR) stars and 
O and B stars. Using far i nfrared (FIR ) and s ub- 
millimcter continuum data, IZvlka et al. es- 
timated that the central cavity (radius ~ 1 pc ) 
contains -400 M Q of dust at T km -40 K, 4 M Q at 
-170 K and -0.01 M Q at -400 K. The luminosi- 
ties emitted by these three phases are — 2 x 10 6 
L(7), 4 x 10 6 L© . and 5 x 10 5 L Q respectively. 
Latvakoski et al.l ( 19991 ) studied the region from 



the Kuiper Airborne Observatory (KAO) in the 30 
/im band with a resolution of 8.5" , and reported 
that structures in the region have temperatures 
ranging from about 70 K to 100 K. They found 
the total luminosity of the components in the in- 
ner region (from their Table 1) to be — 6xl0 6 L Q 
and its mass —2000 M©. 

Surrounding the central cavity, extending from 
1 pc to — 5 pc, is the circumnuclear disk (CND). 
The CND has a mass of several thousand solar 
masses, and rotates with a velocity — 100 Km-s -1 . 
The neutral material in this ring is clumpy, 
dense (n(H 2 ) — 1 3 ~ 7 c m" 3 ), and fragmented 



(Genzel et al. 1982 



■TO % IZvlka et al 



1987) having a filling facto r 



19951 ): iMcOarv et alJ lj200lh 



Wright et al.l (|200lh . There is a clumpy and fila- 
mentary structure known as the mini-spiral that 
connects the CND with the central ionised cav- 
ity, allowing ultraviolet radiation from the central 
ionised cavity to penetrate int o the CND, heat 



ing a nd photoionising the gas (lYusef-Zadeh et al 



2001). 



The kinetic temperatures of the gas in the CND 
and central cavity ranges as high as 200 K - 
300 K in CO(J = 7 - 6) (jBradford et alJ 1200.4 



Martin et all 120041) and -400 K in NH3 (6,6) 



(Montero-Castaho et al.l l2009t Herrnstein et al 



2005). The source of heating remains unclear. 



One explanation suggested is that the ionised cen- 
tral cavity is excited by the radiation from a cen- 
tral cluster of O and B stars, and that at the cav- 



ity edge the surrounding material is likely to be 
shock compressed by the expanding gas around 
the PDRs. Other possibilities include C-shock 
heating, perhaps by magnetic v iscous heating. 



or pos sibly turbulent dissipation. iBradford et al 



([20051 ) suggest MHD (magneto hydrodynamics) 
shock heating with a shock velocity of order 20 
- 30 Km-s" 1 , and a magnetic field 0.3 - 0.5 mG. 
If so, the material in the CND will dissipate its 
orbital energy within a few revolutions, and they 
suggest that its apparent longevity may result 
from refueling by infalling material from the sur- 
rounding clouds. Several authors suggest that the 
ionisation of the molecular clouds in the Galactic 
Centre is d ue to a large flux of cosmic rays heating 



the region dGusten et al.l 



1993 



Yusef-Zadeh et all 



1981 



2007) 



Huttemeister et al 



The CND is surrounded by a number of dense 
molecular clouds and filamentary structures (see 
Figure [T|). At a distance of —39 pc to the 
North East of SgrA* is the Radio Arc consist- 
ing of thermal and nonthermal structures aligned 
almost perpendicularly to the Galactic plane 



(jYusef-Zadeh et al.lll984l ). The Radio Arc seems 



to connect to SgrA* via two concentrations of 
ionised and molecular gas: the Sickle and the 
Arched Filaments. Two clusters of WR and O- 
type stars, the Quintuplet and the Arches Cluster, 
are responsible for the ionisation of the Sickle and 
the Arched filaments, respectively. The Arches 
cluster is slightly in the foreground, about 20 
pc from the Arched filaments, and irradiates 
the several filam ents approximately uniformly 
(|Lang et al.ll200ll ). 

The MSX image at 21.3 /im and the PACS im- 
age at 70 /im trace a shell surrounding the Arches 
and Quintuplet clusters. The shell goes through 
the Arched Filaments, the Lima bean and passes 
through the south of the GO. 11-0. 11 molecular 
cloud f orming a ring of — 44 pc in diametre (Fig- 
ure iBallv et al.l (2 Old) suggest that the radia- 
tion and the stellar winds from the Arches and the 
Quintuplet clusters and the black hole in SgrA* 
may be responsible for the formation of this shell, 
whic h can clearly be seen at t he Spitzer 24 /im im- 
age ( Yusef-Zadeh et al. 20091). and which has also 



been reported bv lMoneti et al] (]2001[ ). 

The objective of this paper is to further con- 
strain the physical conditions in the central 2 par- 
sees of the Galactic Centre, and to help clarify the 



2 



47:00.0 30.0 17:46:00.0 30.0 45:00.0 

Right ascension 

Fig. 1— Composite image of the SgrA* region; MSX-E at 21.3 ^m (blue), 70 fim (green) and SPIRE 350 
\xm (red). The CND appears more prominent at 21.3 fim and 70 /jm due to the lower angular resolution 
of the 350 /im observations. The black cross shows the position of Sgr A* at RA (2000)= 17h 45m 39s. 77 
Dec (2000)= -29° 00' 43.0". A and Q mark the position of the Arches and the Quintuplet clusters of O 
and WR stars, which provide enough photons to ionise and heat the Sickle and the Arched filaments. The 
position of the main molecular clouds (GO. 11-0. 11, M-0.02-0.07 and M0. 25+0. 01) across the Galactic Centre 
are also indicated. Each tic on the Declination axis is an arcminute, corresponding to 2.4 pc. Black contours 
represent fluxes 0.1, 0.45 1.0, 1.97, 2.89, 4.0, 5.0 and 14.0 Jy on the PACS 70 image. 

excitation conditions. We use archival Herschel region to determine the dust properties more ac- 

PACS and SPIRE photometric observations of the curately; many earlier data sets at shorter wave- 
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lengths ( Latvakoski et al. 19991) would have been 
insensitive to a cold dust component. We also in- 
corporate previously unpublished infrared images 
from the Infrared Space Observatory Long Wave- 
length Spectrometer (ISO-LWS), operating in par- 
allel mode; its grating sampling allows us to obtain 
a more precise analysis of the character and struc- 
ture of dust in the central region and the CND. 

2. Observations 



Herschel PACS (jPoglitsch et alj 2001) d ata at 



70 M m and 160 ,um, and SPIRE ([Griffin et al 
20061 ) data at 250 /mi, 350 /mi and 500 /mi, were 
used to measure the dust temperature and the col- 
umn density distribution across the Galactic Cen- 
tre region. The basic data were obtained from 
the public Herschel archive, and were reprocessed 
using HIPE, the Herschel Interactive Processing 
Environment. The data were collected in parallel 
mode and are still rather preliminary in absolute 
accuracy, suffering from calibration errors. PACS 
fluxes off the bright galactic plane are negative, 
indicating that the pipeline corrections are inac- 
curate; moreover there are indications that in this 
extremely bright region some non-linear detector 
responses are uncorrected. We do not apply any 
across-the-board offsets to get positive flux values. 
As a result, the derivative images we discuss be- 
low are reliable over most but not all of the area, 
but the SEDs we fit to the warm central region 
(see Section 4) are accurate. Images are made 
in Jy/pix at 2.0"/pixel and 3.0"/pixel for PACS 
70 /im and 160 /mi, respectively, and in Jy/beam 
for SPIRE data at 250, 350 and 500 /im, with 
beam areas respectively of 395. 0" 2 , 740. 0" 2 and 
1571.0" 2 . Figure Q] covers an area of 28.46 x 29.4 
arcmin 2 centred at RA(2000)= 17h 45m 39.77, 
Dec(2000)= -29° 00' 43.0". For a distance of 8.5 
Kpc, the covered area is ^70.4 x 72.7 pc 2 . 





PACS 




SPIRE 






70.0 


160.0 


250.0 


350.0 


500.0 


pixel (") 


2.0 


3.0 


6.0 


10.0 


14.0 


FWHM (") 


5.9 


11.6 


18.5 


25.3 


36.9 


FWHM (pc) 


0.24 


0.48 


0.76 


1.04 


1.52 



Table 1: PACS and SPIRE pixel and beam sizes 



The image is color-stretched to show more of 
the faint or cool structures, which are labeled. The 
most obvious places where the photometry issues 
make analyses difficult are in the coldest, dark 
regions, for example in GO. 11-0. 11. We white- 
out the corresponding area in the derivative col- 
umn density and spectral index images because 
of the uncertain correction to the negative fluxes 
reported by the current Herschel pipeline. 

The MSX (Midcourse Space Experiment) sur- 
veyed the Galactic Plane at 6 bands from 4.2 to 
26 m ii at a spatial resolution of 18" ( Egan et al.l 
19991 ). The MSX data can also be used to study 
the morphology, dynamics and physics of the ISM 
providing the thermal emissio n from dust grains 
at mid infrared wavelengths (Cohen 19991) . We 
used the MSX Band E data at 21.3 /mi to trace 
the warm dust across the Galactic Centre. 

We also report in this paper previously un- 
published observations of Sgr A*, taken by the 
Long Wavelength Spectrometer (spectral range 
46-180 /mi) on the Infrared Space Observatory 
satellite in parallel mode. ISO was equipped 
with two spectrometers: SWS and LWS, a cam- 
era: ISOCAM, and an imaging photopolarimctcr: 
ISOPHOT to cover wavelengths from 2.5 to 240 
/tin with spatial resolutions of 1.5", at the short- 
est wavelengths, and 85" -90" at the longer wave- 
lengths. Two instruments, LWS and ISOCAM, 
were able to operate simultaneously in parallel 
mode. More than 100,000 individual pointings 
were made in parallel mode in over 17,000 individ- 
ual observations with a total sky coverage of about 
1%. Key to our study, the LWS operate d with ten 



wave lengths in the range 46-180 /mi (jLim et al 
1999), each with bandwidths of about 0.3 /im for 
the five short wavelength detectors (SW) and 0.6 
/tm for the five long wavelength detectors (LW). 
Parallel mode maps are generated by combin- 
ing different raster scans. Several data reduc- 
tion tools were developed in the IDL program- 
ming language and these form the LWS parallel 
interactive analysis package that can be found at: 



|http:/ /jackal. bnsc.rl.ac.uk/ iso /lws / software / soft war e, html | 

Now that the Herschel observations have ob- 
tained longer wavelength coverage, out to 500 /im, 
the earlier ISO data can reliably be combined with 
them to prepare a more comprehensive picture. In 
particular, the spectral information obtained with 
ISO-LWS allows us to determine the peak wave- 
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length and flux of the dust emission with much 
greater precision than is possible with broad-band 
photometry; this in turn allows us to specify the 
conditions of excitation much more accurately. 

3. Dust grain properties; temperature and 
optical depth maps 

Herschel PACS and SPIRE photometry have a 
much higher spatial resolution in the FIR than 
any previous mission. When analysed, the maps 
provide a much improved spatial distribution of 
the fluxes, temperature, spectral index, and H 2 
column density across the Galactic Centre region. 
The MSX 21.3 /jm images provide an important 
measure of the warmer material in the region. Fig- 
ure [T] shows the Galactic Centre as seen by com- 
bined images of MSX-E 21.3 /im (blue), PACS 70 
/im (green) and SPIRE 350 /iin (red). The arched 
filaments, the Radio Arc and Sgr A* are clearly 
visible at 21.3 /im and 70 /im, while SPIRE 350 /im 
traces the very cold and dense molecular clouds. 
Herschel resolves the central cavity and the sur- 
rounding CND at 70, 160 and 250 /im; the lower 
SPIRE 350 and 500 /im resolutions (FWHM= 
25.3" and 36.9" respectively) are not quite able to 
resolve the cavity. 

The observed continuum emission from the 
Galactic Centre region peaks at about 66 /im, in 
the range of ISO-LWS detector SW3 (Figure HJ. 
As shown in Figure [2j the shape of the spectral 
energy distribution (SED) cannot be fit by a sin- 
gle grey body curve. In the ensuing discussion, we 
present our procedure for decomposing the shape, 
which depends on the temperature, the number of 
grains along the line of sight, and their emissivity 
and spectral index, (3. At the longest wavelengths, 
the grain properties determine the shape of the 
SED. We assume a grain size of a ~ 0.1 /im. The 
temperature and the spectral index, /3, distribu- 
tion across the Galactic Centre was then obtained 
by fitting the PACS and SPIRE wavelengths with 
a grey body curve: 



OxB, (T) x (1 - e - T " 



(1) 



where f2 is the solid angle, T is the temperature 
and t v is the optical depth for a frequency v, given 
as: 



To is the bibitem optical depth at uq, the bibitem 
frequency, which was chosen to be c/100 /jm. The 
shape of the curve is only mil dly sensitive to the 



choice of v (jLewis et al. I boooh . 



The SED for each pixel in the Herschel image 
was fitted after smoothing the images and adjust- 
ing the pixel sizes uniformly to the resolution of 
the SPIRE 500 /im channel band. The H 2 column 
density was then derived assuming a gas-to-dust- 
ratio of 100. We did not include the PACS 160 /on 
image data with the other four Herschel bands in 
the calculation of the temperature and the column 
density maps. 

At the brightest regions, the flux values from 
this detector did not fit any physical SED, per- 
haps due to a poor correction of low level stripes 
in the data and / or errors on the background sub- 
traction. Moreover, in the low intensity regions 
the detector map has negative valu es possibly due 



to un corrected gain adjustments ((Bernard et al 



TO 



(2) 



2010). Its omission will not affect the conclusions 
significantly; subsequent pipeline reductions will 
hopefully address these issue. 

Figures |31 |U and [5] show the resultant tem- 
perature, molecular hydrogen column density and 
/3 maps, respectively; north is up in all cases. 
The white circle on the temperature map, cen- 
tred at Sgr A*, shows the average size of the ISO- 
LWS beam which covers, approximately, the cen- 
tral cavity and the inner 1 pc of the CND, for a 
distance to the region of 8.5 Kpc. Figure [7] shows 
the PACS 160 /im map; although the absolute 
calibration is defective, as noted, its higher res- 
olution provides an improved picture of the CND 
than does either the temperature or column den- 
sity maps, which have been reduced to the SPIRE 
500 /xm beam size and are therefore unable to re- 
solve the CND and cavity structure. 

3.1. The temperature map 

The temperature map (Figure [3]) shows the 
most distinguishing details of the various struc- 
tures in this complex region. As expected, due to 
the strong UV radiation and the stellar winds from 
the WR and the OB stellar clusters placed at the 
central cavity, as well as the electron heating due 
to magneto-hydrodynamic turbulences that take 
place in accretion sta tes around the central black 
hole (|Liu et all l201lh . SgrA* shows the highest 
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temperature values , where Tdust > 60 K. 

This map beautifully resolves the Arched 
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Fig. 2.— (Top) Integrated fluxes in the ISO-LWS 
beam for each detector (black triangles) centred 
at the position of Sgr A* . The continuum is fitted 
by the sum of three different dust components: 
(blue) a very hot (Tdust — 90 K) with a very low 
density component (n# 2 ~ 2.8 x 10 3 cm' 3 ) aris- 
ing from the central cavity, and two different dust 
component coming out from the inner part of the 
CND: (green) a hot component (Tdust = 45 K) 
with a density (n# 2 ~ 1.19 x 10 4 cm~ 3 ) and (red) 
a cold component (Tdust — 23.5 K) with very high 
density (n# 2 ~ 2.41 x 10 4 cm -3 ). The grey line 
represents the total contribution. (Bottom) The 
continuum extended to longer wavelengths. The 
black squares represent the Herschel fluxes inte- 
grated on the ISO-LWS beam at each detector: 
PACS at 70 /zm, and SPIRE at 250, 350 and 500 
yum 




TEMPERATURE (K) 

Fig. 3. — Temperature map calculated with the 
PACS 70 fim and SPIRE 250, 350 and 500 /rai. 
The ISO-LWS beam size is shown with a white 
circle centred at the position of Sgr A* . The black 
axis represents longitude (from top to bottom) 
and latitude (from left to right) in Galactic co- 
ordinates. (1) shows the position of the Arched 
filaments, (2) the Sickle, (3) the Radio Arc and 
(4) points at the position of the molecular cloud 
GO. 11-0.11. The size of the image is - 24 x 24 
arcmin 2 (~ 59.4 x 59.4 pc 2 , for a distance of 8.5 
Kpc). 



Filaments and Sickle with dust temperatures 
Tdust ~ 35 K, and the Radio Arc, with signifi- 
cantly lower temperatures than the Arched Fil- 
aments: Tdust ^28 K. As expected, the coldest 
region in the temperature map occurs in the giant 
molecular cloud GO. 11-0. 11, whose temperatures 



arc T 



dust 



< 18 K. As we discuss in more detail 
below and in Figure [SJ these results are in agree- 
ment with the ISO-LWS images, in which the 
Arches (for example) fade away as th e wavelength 



chang es from 46 fxm to 178 fxm. iMorris et al 
(|l995j ) used the KAO at 50 and 90 fxm to map 
this region. They infer a dust temperature for 
the Arches of 45-55 K using a /3 = 1, noticeably 
warmer and inconsistent with our result, in some 
part because (3 is actually closer to 2. 
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Fig. 4. — H2 column density map calculated with 
the PACS 70 urn and SPIRE 250, 350 and 500 
//m. The white square represents the area cov- 
ered by Figure [7] The scale was set to a range of 
(2 - 18)xl0 22 cm" 2 . GO. 11-0.11 molecular cloud 
(4) shows very high H2 column density values 
N(H 2 )~(30 - 80)xl0 22 cm" 2 . The whole molecu- 
lar cloud and the surrounding areas appear like a 
white bulk because the scale was limited to values 
N(H 2 )<18xl0 22 cm~ 2 for abetter definition of the 
lowest column density areas on the map. The size 
of the image is ~ 24 x 24 arcmin 2 (~ 59.4 x 59.4 
pc 2 , for a distance of 8.5 Kpc). 



Colgan et al. I (|l996l) reported that the more 
eastern filament (called E2) might even be higher 
than 70 K. The new Herschel - ISO results show 
that there is considerable knotty structure along 
each of the filaments and elsewhere across the re- 
gion, but that average dust temperatures are more 
accurately closer to 35 K. 

3.2. The column density map 

The column density map (Figure 0]) reveals the 
presence of numerous dense fea tures most of which 
have been previously reported. iBallv et al.1 (l2010h 
discuss many of them based on their SHARCII 
0.35mm and BOLOBAM 1.1 mm continuum im- 
ages with a resolution of 9" and 33' respectively. 



Figure [4] shows only about a factor of ten vari- 
ations across most of the region, with few distinct 
structures being defined by their high column den- 
sities. The exception is the well defined, very dark 
cloud GO. 11-0. 11 with very low dust temperature, 
between 12 K and 18 K, and whose peak column 
densities exceed N(H 2 )> 3.0 x 10 23 cm~ 3 . By con- 
trast, the total luminosity varies like T 4+/3 across 
the region, in which j3 changes from about 1 to 3, 
and the dust surface brightness varies by a factor 
of about 300 even in the less dense structures. The 
column density map closely anti-correlates with 
the temperature map, with denser regions being 
colder. 

As a result of the low fluxes, as noted earlier, 
the PACS photometry in the area of GO. 11-0. 11 
has uncorrected negative values, and the density 
and P values are approximately only. In Fig- 
ures |4] and [5j the scale has been intentionally 
stretched to highlight the less prominent struc- 
tures at the expense of saturating the GO. 11-0. 11 
area). The cloud itself is more accurately defined 
by the SPIRE emission (see Figure [TJ in red), and 
indeed the images reveal considerable knots and 
substructure across the source. The southerly ex- 
tension to the source, visible in the three derived 
maps, is an artifact of the negative photometry in 
the low flux region; SPIRE emission clearly shows 
that the low flux is due to an absence of material 
rather than cold high density dust. We can nev- 
ertheless estimate a lower limit to the total mass 
of G0.ll-0.il as lxlO 6 M Q . 

3.3. The spectral index map 

The map for /3 is shown in Figure [SJ it is also 
anti-correlated with the temperature map (Fig- 
ure [3|). SgrA*, for example, the hottest region, 
shows the lowest values of /3~ 1.0-1.3. The Arched 
Filaments and Sickle show values in the range 

1.3-1. 8 and the range of (3 values in the Radio 
Arc is /3~1.7-2.2. The G0.ll-0.il molecular cloud 
shows extremely large values of (3 ^3.5. This gen- 
eral trend is expected; dust grains in cold clouds 
tend to be larger in size and t herefore have larger 
than grains in warm clouds (jButler et aL 2009; 
Hirashita etldl l2007t iDupac et al.l l2003h . The 
strong and rapid spatial variations in the Galactic 
Centre, coupled with the Herschel maps, provide 
one of the first clear looks at how this feature of 
grains manifests itself in a complex region. 
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SPECTRAL INDEX (beta} 

Fig. 5. — Spectral index map calculated with the 
PACS 70 fun and SPIRE 250, 350 and 500 fim. 
(1) shows the position of the Arched filaments, (2) 
the Sickle, (3) the Radio Arc and (4) points at the 
position of the molecular cloud GO. 11-0. 11. The 
size of the image is ~ 24x24 arcmin 2 (~ 59.4x59.4 
pc 2 , for a distance of 8.5 Kpc). 



Since the total luminosity varies roughly like 
T i+ p , this parameter also helps sort out the ex- 
citation and energetics of the various structures. 
The anticorrelation between temperature and den- 
sity applies to most of the other dense regions of 
the maps as well, which, being dominated by dust 
properties, are coldest where they are most opaque 
(Figure UJ; for example through the Arched Fila- 
ments, Sickle, the Radio Arc and SgrA*. 

3.4. The distribution of the very small 
grains (VSG) 

The O and WR stars from the Arches and 
the Quintuplet clusters (Figure [T]) provide enough 
photons to ionise and heat Sic kle and the Arched 
filaments (jCotera et al.l Il996l ) and also generate 
strong winds that could contribute to the deple- 
tion of the largest dust grains whose population is 
traced by the N(H2) map: The largest dust grains 
can be destroyed due to the strong winds and the 




MSX-E (xio 4 W/m2/sr) 

Fig. 6.— MSX-E image at 21.3 fan. (1) shows the 
position of the Arched filaments, (2) the Sickle, 
(3) the Radio Arc and (4) points out the position 
of the molecular cloud GO. 11-0. 11. The size of the 
image is ~ 24 x 24 arcmin 2 (~ 59.4 x 59.4 pc 2 , for 
a distance of 8.5 Kpc). 



UV radiation from the massive O and WR stars 
from the Arches and Quintuplet clusters. 

The MSX-E map at 21.3 fim (Figure [BJ can 
be used to trace the population of very small 
grains (VSG) throughout the Galactic Centre re- 
gion (this band is more reliable than the features 
observed at shorter infrared wavelengths which 
can be affected by extinction). VSGs are small 
(with sizes of ~ 0.01 p,m), and a bsorption of 
one UV photon ( Pomares et al.|[2009h can increase 
their temperature up to 80 K, prompting them to 
re-emit at these shorter wavelengths (A ~ 24 fim) 



([Desert et al.l 1990; B oulanger et all 119851) . The 



MSX-E map correlates well with the temperature 
map, showing that the population of the small 
grains is distributed through the regions where the 
population of the largest grains has been depleted, 
i.e. the Sickle and the Arched filaments. This 
effect can also be observed in SgrA*: PACS 160 
fim, Figure [7], shows very low emission at this far- 
infrared wavelength in the central cavity of SgrA* , 
but which is surrounded by bright emission from 
the CND. The emission in the inner parsec, as 



S 



Flux (Jy) 5.8 



9.3 13 16 20 23 




Right ascension 

Fig. 7. — SgrA* image from PACS 160 /im. The 
map resolves the CND at 160 /jm (black con- 
tours) surrounding the central cavity, which is 
well-traced by the MSX-E contours (white) at 21.3 
/im. The MSX-E intensity peaks in the centre 
of the cavity at 3.89xl0~ 3 W/cm 2 /sr decreasing 
outwards to 1.0xl0~ 3 W/cm 2 /sr. The size of the 
map is ~ 3.3 x 3.4 arcmin 2 (^ 8.2 x 8.4 pc 2 , for a 
distance of 8.5 Kpc). 



traced by MSX-E at 21.3 /im, appears confined to 
the central cavity where the O and B stars from 
the central cluster produce a strong radiation field, 
and where stellar winds heat up the small grains 
and deplete large grains. Within the central 2 par- 
sees, the largest grains are found remaining in the 
CND itself, where they are shielded from the cen- 
tral strong UV radiation; their emission is traced 
by Herschel PACS+SPIRE at far infrared wave- 
lengths A > 60 /im. 

4. The Temperature Components of the 
CND: Fitting the ISO LWS continuum 

ISO-LWS obtained images of the central region 
of the galaxy operating in its parallel mode, one 
image simultaneously from each of the ten ISO- 
LWS detectors, with the spectrometer's grating set 
at a fixed position near the centre of its scan range. 



The ISO LWS parallel mode maps provide us with 
complementary wavelength coverage just where it 
is most needed, at the peak of the dust emission, 
in the range of 46 to 180 /im that can be used 
to defined the SED of SgrA* at far infrared wave- 
lengths. Unfortunately, with a beam of ~ 85" and 
a pixel size of 42" , these ISO maps are unable to 
resolve as many structures as do higher resolution 
Herschel PACS and SPIRE data. 

Figure [5] shows the resultant Galactic Centre 
images in each of the ten bands. The images 
are centred at RA (2000)= 17h 45m 56s.38, Dec 
(2000)= -28° 55' 15.24" an cover and area of ~ 
22.8x29.6 arcmin 2 (~ 59x77 pc 2 ). In order to 
obtain accurate photometry from these images, 
an extended source flux correction factor is re- 
quired, as described in the ISO Users Manual. For- 
tunately we also had pointed observations with 
which to confirm the correction procedure (see 
(jWhite et al.l 120111 ). for a description of the ISO 
pointed photometric and spectroscopic observa- 
tions). The observed flux in Jy is: 



S v [Jy] = F v x //(fi x 10 6 )[MJy ■ sr' 



(3) 



and the extended source correction factors (/) 
and t he solid angle O are given in Table 2. Be- 
sides, lCasassus_eLa]J (|2008l ) reported that the ISO 
LW1 at 102.42 /im intensities were 27% higher 
than IRAS 100 /im for intensities bellow 1000 
MJy-sr -1 and 40% higher for intensities above 
1000 MJy-sr" 1 . This was in agreement with 



Chan et al.l (|2001l ). This offset on the ISO par- 
allel mode data was calculated by comparing the 
fluxes at different position along the Galactic Cen- 
tre observed in primary mode with the fluxes ob- 
served in parallel at the same position for the 
ten ISO LWS detectors. The correction factors, 
Q{Flux para u e i/Flux primar y), are also given on 
Table 2. 

Despite the correction factors were applied, the 
LW2 fluxes seem to be still offset in our data by 
10% above the continuum level. Due to this prob- 
lem, the LW2 flux was omitted when fitting the 
ISO-LWS SED in Figure 0. 

SgrA* is very bright at the shortest wave- 
lengths, and clearly seen in Figure [5] The SW 
detectors (46.2 - 84.79 /im) also show the north- 
ern part of the Radio Arc. Detectors SW4 75.69 
/im, SW5 84.79 /im, LW1 102.42 /im and LW2 
122.2 /im show SgrA* as a single object about ^4 
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Fig. 8.— Galactic Centre observed at each ISO-LWS parallel mode wavelength: 46.222 /an (SW1), 56.2033 
/im (SW2), 66.1173 /an (SW3), 75.6989 /im (SW4), 84.7977 /im (SW5), 102.425 /im (LW1), 122.218 fim 
(LW2), 141.809 /im (LW3), 160.554 /im (LW4) and 177.971 fj,m (LW5). Black contours on SW3 are PACS 
70 /im convolved with the SPIRE 500 /im beam size (36.9") and the black contours on LW4 correspond to 
PACS 160 /im emission with a beam size of 36.9" . The black circle on the SW1 map represents the ISO-LWS 
beam (~ 3.5 pc) centred at the position of SgrA*. 



Detector 


X(fim) 


/ 


Q x 10 6 [sr] 


Q 


SW1 


46.22 


0.8704 


0.1140 


1.082 


SW2 


56.20 


0.8677 


0.1321 


0.899 


SW3 


66.11 


0.8421 


0.1410 


0.889 


SW4 


75.69 


0.7334 


0.1223 


0.823 


SW5 


84.79 


0.6878 


0.1163 


0.878 


LW1 


102.42 


0.6758 


0.1111 


0.930 


LW2 


122.21 


0.6734 


0.1128 


0.823 


LW3 


141.80 


0.6035 


0.0935 


0.948 


LW4 


160.55 


0.5411 


0.0904 


0.992 


LW5 


177.97 


0.4596 


0.0833 


1.004 



Table 2: Table of extended source correction (/) 
and effective solid angle of the beam (f2) for the 
different LWS detectors. 



pc in radius that includes SgrA* and extends to M- 
0.02-0.07 molecular cloud (Figure [lj in the south. 
At 160.55 fj,m (LW4) and 177.97 /im (LW5), the 
most conspicuous feature is the cold and dense 



molecular cloud M-0.02-0.07 (see also Figure [lj. 

The spectral energy distribution of Sgr A* 
within the ~85" ISO-LWS beam is shown in Fig- 
ure [3J Each data point on Figure [2] represents the 
integrated flux on the ISO-LWS beam centred at 
the position of SgrA* at each wavelength. We note 
that at a distance to the region of about 8.5 Kpc, 
the ISO-LWS beam size of 3.5 pc includes the cen- 
tral cavity and the inner ~0.8 pc of the CND, not 
quite large enough to completely include the full 
CND observed in the Herschel maps. In order to 
model the ISO-LWS continuum spectrum across 
the wavelength range ~40-200 /im, we assumed 
typical dust properties: dust grains of 0.1 /im in 
radius and a bulk density of 2.5 (g-cm~ 3 ). These 
values are characteristic of silicate gr ains that con- 



tain H 2O inside their structure (e.g. iPollack et al 
([l99j)Pollack et al. 1994). The dust opacity is 
defined as: 

k\ ■ p = (4) 



where KA(cm 2 -g x ) is the dust opacity, p(g-cm 3 ) 
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is the bulk density of the grains and K cons t is a 
proportionality constant. For A =100 /urn and 
K CO nst= 0.8 the dust opacity in the central cav- 
ity is kioomwi ~ 92 cm^g" 1 , characteristic of dust 
grains with thin ice mantles, while the opacities 
at the CND are -160 and 200 cm 2 /g for the hot 
and the cold component, respectively. These val- 
ues are typical of dust grains with thick ice layers 



(|Ossenkopf fc Henning Ill994l ) 



The best model that we found able to reproduce 
the LWS spectral continuum consists of three- 
temperature components: an SED arising from 
the central cavity of r=l pc in radius, and a com- 
bination of a warm plus cold SED arising from 
the inner 0.86 pc of the CND. Figure 0] shows the 
steep gradient in column density at the position of 
SgrA*, from values of NH 2 — 3 x 10 22 cm" 2 at the 
centre of the cavity to higher values ~ 1 x 10 23 
cm~ 2 at the CND. The temperature map (Fig- 
ure [3]) likewise shows the largest values just on 
the centre , T, 



dust 



> 



60 K, surrounded by colder 
dust. The SED from the cavity can be fitted by a 
grey-body curve defined by a dust temperature of 
Td us t= 90 K, /3=1.23 and a low column density of 
N(H 2 )= 8.73 x 10 21 cm" 2 . 

As noted, the SED of the CND requires two 
components in order to fit the whole ISO-LWS 
continuum: a hot component of Tdust— 44.5 K, 
/3=1 .35 and a column density of N(H 2 )= 3.18 x 
10 22 cm -2 . Significantly, a cold component also 
has to be introduced in order to fit the longest 
wavelengths of the ISO-LWS data. The cold com- 
ponent can be fitted by a r Tdust= 23.5 K modified 
blackbody with /3=1.40 and a higher column den- 
sity N(H 2 )= 6.39 x 10 22 cm -2 . This massive new 
cold component, previously unreported in obser- 
vations done only at shorter wavelengths, is re- 
confirmed by the longer wavelength Herschel pho- 
tometry. This simple, three component model re- 
produces very well the observed continuum emis- 
sion at the ISO-LWS spectrum (Figure [5]) and the 
SPIRE and PACS 70 ^m data. The luminosity 



2.2 x 10 b L Q , 
- 1.0 x 10 5 L @ , 
~ 1.7 x 10 3 M Q , 
m cnd ^ 1-56 x 10 4 Mq and Mg$ D - 3.15 x 10 4 



of each component is 
^CND ~ 1-4 x 10 6 L Q 
and the total gas mass 

hot _ , i cfi „ 1 r>4 



^cavity 

ana i^ CND 



M 



cavity 



Mq . Table [3] lists the parameters from the fit to 
the ISO-LWS continuum. 

Th e flux errors on t he ISO-LWS data are 10%- 
15% (jGrv et al.l [2003) . Assuming an error of 





Cavity 


CNDft ot 


CND coid 


radius(pc) 


1.0 


0.86 


0.86 


P 


1.23 


1.35 


1.40 


Kconst. 


0.8 


0.8 


0.8 


KlOO/jm 


92.28 


160.38 


200. 


Tdust{K) 


90.0 


44.5 


23.5 


Pdust {cm' 3 ) 


9.45X10 -23 


4.0xl0~ 22 


8.05xl0~ 22 


TloOjim 


5xl0" 3 


3xl0~ 2 


8xl0 -2 


Nh 2 (cm- 2 ) 


8.73X10 21 


3.18X10 22 


6.39X10 22 


riH 2 (cm -3 ) 


2.8xl0 3 


1.19X10 4 


2.41 xlO 4 


Xdust 


0.01 


0.01 


0.01 


£(£©) 


2.2 x 10 6 


1.4 x 10 6 


1.0 x 10 5 


M (Mq) 


1.7 x 10 3 


1.56 x 10 4 


3.15 x 10 4 



Table 3: Parameters used in order to fit the ISO- 
LWS continuum. 



~12% for each detector, the goodness of the fit 
on Figure [5] is x 2 = 0.89 (LW2 flux was not taken 
into account on the estimation of x 2 )- The CND 
cold component has to be included in the model in 
order to fit accurately the SPIRE fluxes at wave- 
lengths > 250 /im. Models with just one CND 
component can also reproduce the ISO-LWS con- 
tinuum very precisely. However, the model must 
satisfy another constraint set by the fluxes of the 
atomic and molecular lines detected by ISO-LWS 
in this region. We have used a radiative trans- 
fer code as described in Gonzalez- Alfonso et al 



(|2002f ), to simultaneously fit the continuum and 
the spectral lines, obtaining the best results for 
the case of two CND components. We present the 
discussion of the gas properties thro ugh SgrA* in 
a separate paper (jWhite et al.ll201ll ). 



5. Conclusions 

We present new Herschel PACS and SPIRE 
maps of the Galactic Centre, together with pre- 
viously unpublished ISO-LWS 10-band parallel 
mode maps that enable us to obtain an accurate 
calibration for the central parsecs and a detailed 
spectral shape around the wavelength peak of the 
continuum emission. Photometry from the data 
set allows us to model dust in the central 2 pc of 
the Galactic Centre. Combined with MSX-E 21.3 
jum images, we use these new results to trace the 
dust properties throughout the Galactic Centre re- 
gion. 

The ISO LWS continuum in the direction of 
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Sgr A*, extended with Herschel data out to 500 
/jm, is best fit by a model of a central cavity of 1 
pc radius emitting like a greybody with T = 90 K 
and P= 1.23 and with low density n(H 2 )= 2.8 xlO 3 
cm~ 3 . The region is surrounded by the circum- 
nuclear disk represented by two different compo- 
nents: a hot component with a dust temperature 
of 44.5 K and n(H 2 )= 1.19xl0 4 cm" 3 with /3= 
1.35, and a cold component with T dust — 23.5 K, 
n(H 2 )= 2.41 xlO 4 cnr 3 and /3= 1.4. The two com- 
ponents indicate that the CND is not a single uni- 
form body, but (as suggested by the images and 
previous authors) has numerous clumps, gaps, and 
substructures but whose average properties mani- 
fest these two basic behaviours. This simple model 
is able to reproduce the observed ISO-LWS con- 
tinuum as well as the SPIRE data at 250 /xm, 350 
/xm and 500 /im. The MSX-E map traces the emis- 
sion of the very small grains population through 
the Radio Arc, the Sickle, the Arched Filaments 
and confined in the central cavity in Sgr A* . These 
spatial structures are now directly detected in the 
far infrared, and, with the full SED available, their 
dust features are all characterized. As expected, 
the largest grains are distributed throughout the 
CND and correspond to the densest and coolest 
molecular clouds. Unresolved issues with the Her- 
schel photometry in the current dataset leaves the 
absolute accuracy of the coldest, densest regions 
(like GO. 11-0. 11) approximate. 

Dust opacity properties in the central cavity are 
typical of grains with a thin ice layer; while the 
dust properties in the CND, with lower tempera- 
ture and highest densities, is more typical of dust 
grains with thick ice mantles. 

The total integrated continuum flux around the 
inner edge of the CND and the central cavity pro- 
vides a total mass of M~4.9 x 10 4 M©. The lumi- 
nosities obtained are Jj ca vity ~ 2.2 x 10 6 L© and 
Lcnd ~ 1-5 x 10 6 Lq in the central ^2.0 pc in 
Sgr A*. 
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made. We are very grateful to Tanya Lim for 
her suggestions on the treatment of the ISO LWS 
parallel mode data. This work was supported in 
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